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Structure formation in the universe can produce high energy gamma-rays from shock-accelerated 
electrons, and this process may be the origin of the extragalactic gamma-ray background (EGRB) 
as well as a part of the unidentified sources detected by EGRET in the GeV band, if about 5% of 
the kinetic energy of the shock is going into electron acceleration. However, we point out that the 
production of gamma-rays may be severely suppressed if the collapsing matter has been preheated by 
external entropy sources at the time of gravitational collapse, as can be inferred from the luminosity- 
temperature (LT) relation of galaxy clusters and groups. We also make a rough estimate of this 
effect by a simple model, showing that the EGRB flux may be suppressed by a factor of about 30. 
Hence structure formation is difficult to be the dominant origin of EGRB if preheating is actually 
responsible for the observed anomary in the LT relation. The detectable number of gamma-ray 
clusters is also reduced, but about 5-10 forming clusters should still be detectable by EGRET all 
sky, and this number is similar to that of the steady and high-latitude unidentified sources in the 
EGRET catalog. The future GLAST mission should detect 10 2 -10 3 gamma-ray clusters of galaxies 
even if the intergalactic medium has been preheated. 



I. INTRODUCTION 



It is widely believed that the observed structures in 
the universe have been produced via gravitational insta- 
bility. Currently the most successful theory of structure 
formation is the cold dark matter (CDM) scenario, in 
which the structures grow hierarchically from small ob- 
jects into larger ones. When an object collapses gravita- 
tionally and virializes, the baryonic matter in the object 
is heated by shock waves up to the virial temperature, 
and particles are expected to be accelerated to high en- 
ergy by shock acceleration. High energy gamma-rays are 
then expected to be produced during structure forma- 
tion, via inverse-Compton scattering of cosmic microwave 
background photons by high energy electrons. Recently, 
Loeb & Waxman |Q speculated that this process can be 
the origin of the extragalactic gamma-ray background 
(EGRB) observed at - 1-100 GeV [§, if about 5% of 
the post-shock thermal energy is going into electron ac- 
celeration. Totani & Kitayama (Paper I ||) has shown 
that, if this is the few tens of forming clusters 

should have already been detected by the EGRET ex- 
periment which has performed an all sky survey in the 
GeV band f§, and a part of the unidentified EGRET 
sources can be accounted for (see also Waxman & Loeb 
§)• 

All these analyses are based solely on hierarchical 



structure formation in the CDM universe, but recent x- 
ray observations of clusters and groups of galaxies have 
shown that the x-ray properties of these objects can- 
not be explained by the above simple picture of hier- 
archical structure formation alone. It is well known 
that the luminosity-temperature (LT) relation of clus- 
ters and groups is considerably different from what is 
expected from the self-similar model predicted by hierar- 
chical structure formation (e.g., ||r||). The most popu- 
lar explanation is that the intergalactic medium has been 
preheated by external entropy sources such as supernovae 
or active galactic nuclei up to a temperature of about 1 
keV (e.g., §||). 

If this is the case, gamma-ray production from struc- 
ture formation should be significantly suppressed, be- 
cause the external entropy impedes the gravitational col- 
lapse and weakens the shock, resulting in decreased shock 
heating and softer spectra of accelerated particles. It 
should be noted that the objects on which preheating 
has the most significant effect (M ~ 10 14 M Q and T ~ 1 
keV) are those which are expected to produce most of the 
EGRB photons. In this letter we try to make a quantita- 
tive estimate of the preheating effect on the production 
of gamma-rays from structure formation. 

Throughout this paper, we assume a CDM universe 
with the density parameter f2o = 0.3, the cosmolog- 
ical constant J7a = 0.7, the Hubble constant h = 



1 



Ho/(100km/s/Mpc) = 0.7, the baryon density param- 
eter J7b = 0.015/i~ 2 , and the density fluctuation ampli- 
tude as — 1. These parameters are favored from various 
recent cosmological observations. 



II. EFFECT OF PREHEATING ON STRUCTURE 
FORMATION 

The formulation for calculating the EGRB flux and 
source counts of gamma-ray clusters when there is no pre- 
heating has been given in Paper I. Here we describe the 
modification to include the effect of preheating. There 
are two important effects of preheating on the production 
of high energy gamma-rays from structure formation. 
The first is that the external entropy leads to the virial- 
ization of collapsing gas at a larger radius with smaller 
infalling velocity and smaller kinetic energy of shock com- 
pared with the no-prcheating case. The second is that 
the external entropy results in a smaller Mach number 
and a shock which is no longer ideally strong, and hence 
the energy spectrum of accelerated electrons should be 
significantly softer than that without preheating. 

It is not easy to predict realistic density profiles for 
preheated, collapsing gas without recourse to detailed 
numerical simulations, even in the spherically symmet- 
ric, one-dimensional case (e.g., [ [l2||Tl| l ) . In the following 
discussion, we do not inquire about the actual distribu- 
tion of the baryonic matter within the dark matter halo, 
and the relevant physical quantities are to be interpreted 
in a volume-averaged sense. It should be noted, how- 
ever, that the production of gamma-rays is less sensitive 
to the gas density profile than, e.g., the x-ray luminosity 
of thermal bremsstrahlung because the target photons 
for inverse-Compton scattering into gamma-rays are the 
CMB photons whose density is universal. Thus we be- 
lieve that the following simplified treatment is an ade- 
quate first approximation for the estimate of gamma-ray 
production, although it may be too simple to calculate 
accurately the x-ray luminosity and temperature of clus- 
ters of galaxies. 

We calculate the virial radius (rh) and the velocity of 
collapsing, preheated baryonic gas (Vh) at virialization 
by a simple model naturally extended from the standard 
spherical collapse model (e.g., |l3|). Let us start from the 
energy conservation equation of baryonic gas with mass 
Mb = (S!b/Oo)M originally embedded in a dark halo 
with mass M: 
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where r is the characteristic radius of the collapsing gas, 
and V — dr/dt is the infalling velocity. We ignore the 
heating or cooling during the collapse and assume the col- 
lapse proceeds adiabatically, and then the effect of pre- 
heating is to add the internal energy MbO 2 /[T(T — 1)], 
where a = (dP/dps) 1 ^ 2 is the sound velocity. The pres- 
sure P is related to the density through the entropy pa- 
rameter K as P = Kp T B , and the baryon gas density is 
Pb = Ms/(47rr 3 /3). The right hand side of this equa- 
tion is the total energy of this system, estimated at max- 
imum expansion (r = 2r v j r ), where r v i r is the virial ra- 
dius when there is no preheating. The sound velocity 
a — (dP/dps) 1 / 2 is related to that at maximum expan- 
sion a a as a = (r/2r vir ) _1 - 5 ( r ~ 1 )a . 

The virialization of the system is expected to occur at 
a radius of which is larger than r v ; r because of pre- 
heating. We estimate this by simply extending the virial 
theorem for this system, resulting in the following equa- 
tion at virialization: 
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Therefore, Vh and are determined by solving equations 
|l| and E| For the case of T = 5/3, the result is 
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The above model is valid only so long as rh < 2r v i r ; 
haloes with > 2r v ; r cannot collapse from the point 
of maximum expansion (turn around) and we assume 
that such haloes have radii = 2r vn . and neither shocks 
nor gamma-rays are generated. The original virial radius 
r v i r can be calculated by the standard spherical collapse 
model for a dark halo of mass M collapsing at z, and 
then we can calculate the shock energy as the kinetic en- 
ergy given to infalling baryonic gas, (3/4)MbV 2 , which 
is reduced compared with the no-preheating case. 

Next we calculate the Mach number of the shock and 
spectral index of shock-accelerated particles. In the rest 
frame of the infall (i.e., undisturbed) gas, the propagation 
of the shock can be regarded as a problem of supersonic 
piston moving with the velocity of infall gas, Vh, mea- 
sured in the cluster rest frame. The shock velocity V s 
(i.e., the upstream gas velocity towards the shock in the 
rest frame of the shock front) is given as (e.g., pl[): 
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and the upstream Mach number is given by M. — 
V s /a. The particle index is given by a = 
-d(log7V e )/d(log7 e ) = (r + 2)/(r - 1), where r = {T + 
l)/[(r — 1) + 2/A4 2 ] is the compression ratio and j e is 
the electron Lorenz factor. 

We must determine an appropriate value of the en- 
tropy parameter, K. so that the above model is con- 
sistent with the observed LT relation of clusters and 
groups. We parametrize the entropy parameter as K — 
-^34,0(1 + z) _1 10 34 erg cm 2 g -5 / 3 , assuming the redshift 
dependence of K oc (1 + z) -1 which is consistent with 
the observed LT relation Following Tozzi & Nor- 

man pi]] , we take the parameter .^34,0 ~ 0.8 as a fiducial 
value to be consistent with the LT relation. 

Figure |l| shows the typical parameters of preheated 
haloes as obtained above. It is clear that the objects 
with M ~ 10 14 Mq is seriously affected by preheating, 
with the particle acceleration index much softer than the 
strong-shock limit of a = 2. Therefore cosmological ob- 
jects less massive than ~ 10 14 M© hardly contribute to 
the high energy gamma-ray background in the GeV band. 

As mentioned above, our model may be too simple to 
calculate the x-ray luminosity and temperature of pre- 
heated haloes, since we have ignored the density profile 
within the halo to which the x-ray luminosity is very 
sensitive. In spite of this difficulty, however, our model 
reproduces the observed LT relation fairly well, as shown 
below. A simple scaling relation for the x-ray luminosity 
of thermal bremsstrahlung emission from cluster gas is 
L oc pbT 1 / 2 where T is the gas temperature. It is ex- 
pected that most of the kinetic energy of collapsing mat- 
ter calculated above is eventually converted into thermal 
energy. We estimate the temperature so that the final 
thermal energy is the total of this energy from gravita- 
tional collapse and the external energy of the preheating. 
Then the temperature of the preheated gas is given by 
T h = (V h /V c ) 2 T vil + f im p Kp 1 g- 1 /kB, where V c and T vir 
are respectively the circular velocity and virial tempera- 
ture of a halo without preheating predicted by the spher- 
ical collapse model. Here p is the mean molecular weight 
and ks is the Boltzman constant. The X-ray luminosity 
is then given by L = (r vh ./r /l ) 3 (T /l /T vir ) 1 / 2 i ss , where 
Lss is the x-ray luminosity of the self-similar model for 
non-preheated clusters of galaxies which is a function of 
cluster mass and formation redshift. We used the formula 
of Kitayama & Suto |Q for L$s- 

Then we can calculate the luminosity and tempera- 
ture of a preheated halo with mass M and collapsing at 
redshift zp. Figure § shows the LT relation predicted 
for cosmological objects observed at z b s = 0, compared 
with observations. Objects should have various zf even 
for the same mass, and here we utilize the distribution 



of zp as a function of mass and z b s derived by Lacey 
& Cole |l6j. The thick solid line shows the LT relation 
when the median of the Lacey-Cole zp distribution is 
applied, using the standard value of the entropy param- 
eter: i^34.o = 0.8. The behavior of the model LT curve 
changes abruptly at (L,T) = (10 42 3 erg s" 1 , 0.6 keV), 
and this corresponds to the point at which th becomes 
equal to the maximum expansion radius, 2r v ; r , and hence 
haloes cannot gravitationally collapse. The dashed and 
dotted lines show the dispersion of the LT relation corre- 
sponding to that of zf, encompassing 68% (1 a) and 95% 
(2 a) of the probability distribution of zp, respectively. 
The thin solid lines are for median zf, but with different 
values of if 34,0 = 0.4 and 1.6. Figure || shows that our 
simple model is in reasonable agreement with the data 
when the same entropy parameter as Tozzi & Norman 
is used. 

Now we can calculate the EGRB flux and source counts 
of gamma-ray clusters by the formulation given in Paper 
I, simply replacing the shock energy and particle acceler- 
ation index by those obtained above. In the following we 
will apply the above model to calculate the EGRB flux 
and spectrum, and the expected counts of gamma-ray 
clusters. 

III. EXTRAGALACTIC GAMMA-RAY 
BACKGROUND 

The efficiency of energy injection from the kinetic en- 
ergy of infalling gas into nonthermal electrons by the 
shock acceleration acceleration is parametrized by £ e , 
and we assume 5% injection, i.e., £ e = 0.05 following 
Loeb & Waxman Q and Paper I. It is widely accepted 
that low-energy cosmic rays are accelerated in supernova 
remnants, and its energy injection efficiency is about 1— 
10% for ions. On the other hand, injection into electrons 
is still highly uncertain, both observationally and theo- 
retically. As is well known, the energy flux of cosmic-ray 
electrons observed above the Earth's atmosphere is about 
100 times lower than that of cosmic-ray protons, and it 
may suggest that energy injection is considerably lower 
for electrons than for ions. However, the local cosmic-ray 
flux ratio does not necessarily reflect that at the produc- 
tion site, because of different energy loss and propagation 
processes. Hence, we take £ e = 0.05 as a maximally pos- 
sible value, although it could be rather optimistic. 

Figure || shows our calculation of the EGRB flux and 
spectrum based on the model described above. We have 
used a value of the magnetic field parameter, £b = 10~ 3 , 
which is the ratio of magnetic energy to the total gravita- 
tional energy given to baryonic matter. Magnetic fields 
of £b ~ 10~ 3 arc sometimes observed in intracluster mat- 
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ter, and the EGRB spectrum extends up to ~ 100 GeV 
if £b > 10~ 5 (Paper I). Note that this parameter deter- 
mines only the maximum photon energy of the EGRB 
spectrum, and the EGRB flux is rather insensitive to 
this uncertain parameter. It can be seen that the EGRB 
flux is severely decreased by a factor of about 30 for pre- 
heating of ^34,0 ~ 0.8, which fits best to the observed 
LT relation. Even if we use a relatively small value of 
-^34,0 = 0.4, the EGRB flux is about one order of mag- 
nitude smaller than that observed. Therefore, if the in- 
tergalactic medium is actually preheated, the structure 
formation is very unlikely to be the origin of EGRB. It 
should be noted that we have already assumed a rela- 
tively large efficiency of shock energy injection into elec- 
tron acceleration, £ e = 0.05, and hence we cannot take 
the option of increasing this parameter to save this hy- 
pothesis. 



IV. FORMING GAMMA-RAY CLUSTERS OF 
GALAXIES 



The top panel of Fig. [| shows the calculation of the 
source counts for gamma-rays from structure formation. 
Compared with the case of no preheating, the source 
counts are decreased by a factor of 10 above the EGRET 
sensitivity limit. See the bottom panel of the figure 
for the mean values of physical quantities of gamma- 
ray emitting objects (mass, redshift, and angular radius) 
brighter than a given flux. The mean spectral index for 
gamma-ray clusters detectable above 100 MeV is rather 
insensitive to the flux, increasing as a = 2.20 to 2.26 
from the sensitivity of the EGRET to that of GLAST, 
because only objects with sufficiently hard spectrum can 
be observed at GeV energies. 

For the canonical value of the entropy parameter, 
^34.0 — 0-8, the number of gamma-ray clusters de- 
tectable by EGRET is about 5 in all sky. If we take 
a relatively small value of ^34,0 = 0.4 to account for 
the model uncertainty, the number is increased to 14. 
These numbers are smaller than those reported by Paper 
I because of the preheating effect, but it is interesting 
to note that the number of steady (i.e., unlikely to be 
AGNs) unidentified sources of the EGRET catalog pre- 
sented in Gehrels et al. @ is 7 for |6| > 45° (~ 24 ± 9 
in all sky), which is not very different from our result. A 
significant part of these high-latitude, steady unidenti- 
fied EGRET sources could be explained by dynamically 
forming gamma-ray clusters. 



V. CONCLUSION 

We have shown that the preheating of intergalactic 
medium, which may have occurred as indicated from 
the observed x-ray properties of clusters and groups of 
galaxies, significantly suppresses the production of high- 
energy gamma-rays from structure formation. If pre- 
heating is actually responsible for the steepening in the 
LT relation, structure formation cannot be the domi- 
nant origin of EGRB. The number of discrete sources 
detectable by the EGRET is also decreased by preheat- 
ing, but 5-10 gamma-ray clusters could still be observ- 
able all sky, which may constitute a part of the uniden- 
tified sources. Even if the preheating effect is profound, 
the future GLAST mission may detect about 100-1000 
gamma-ray clusters, and it may be used to probe the pre- 
heating processes in the intergalactic medium as well as 
the dynamical processes of structure formation. 
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FIG. 1. The virial radius (ru) and velocity (Vh), Mach number, and particle acceleration index (a) of gravitationally bound 
objects under the effect of preheating, as a function of dark halo mass. The radius and velocity are given as the ratios to 
the original virial radius (r v ir) and velocity (V c ) in the case of no-preheating. The solid line is for the objects forming at 
zf = 0, while the dashed line for zf = 1. The entropy parameter is assumed to be K — ^34,0(1 + 2) _1 10 34 erg cm 2 g -5 ^ 3 with 
#34,0 = 0.8. 
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FIG. 2. The x-ray luminosity-temperature (LT) relation of galaxy clusters and groups. The thick-solid line is the model 
prediction for objects observed at z b s = 0, using the median of the Lacey & Cole distribution function for the formation 
redshift (zf)- The entropy parameter is assumed to be K = if34,o(l + z)~ 1 10 34 erg cm 2 g~ 5 ' 3 with ^34,0 = 0.8. The solid 
circles are indicating the grids corresponding to the cluster masses of 10 12 , 10 13 , 10 14 , and 10 15 M Q . The dashed and dotted 
lines show the dispersion of the LT relation due to that of zf, in which zf is included by a probability of 68% (1 a) and 95% (2 
a), respectively. The two thin-solid lines are the same as the thick-solid line, but for different values of the entropy parameter, 
-^34,0 = 0.4 and 1.6. The dot-dashed line is the prediction of the self-similar model of x-ray luminosity of clusters without 
preheating, assuming zf = z b s = 0. The data points are from Markevitch (cross, Q), Arnaud & Evrard (open circle, jfj), and 
Helsdon & Ponman (open square, j^]). 
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FIG. 3. The spectrum of the cosmic gamma-ray background (EGRB) in the GeV band expected from structure forma- 
tion. The short-dashed, solid, and long-dashed lines are for the entropy parameters of ^"34,0 = 0.4, 0.8, and 1.6, where 
K = A'34.o(l + z)~ 1 10 34 erg cm 2 g~ 5 ' 3 . (-^34,0 = 0.8 best fits to the observed LT relation of galaxy clusters.) The dot-dashed 
line is the result of Paper I without the effect of preheating. All thick lines take into account the absorption of gamma-rays in 
the intergalactic field using the opacity presented in Totani (2000) JlTj while the thin lines do not. (Reproduction of gamma-rays 
is not taken into account in either case, see Paper I). The observed data are from Sreekumar et al. (1998) 0. 
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FIG. 4. The upper panel: the cumulative flux distribution of gamma-ray emitting clusters of galaxies. The four curves 
are for different values of the entropy parameter, K — i^34,o(l + ^) _1 10 34 erg cm 2 g _5/ ^ 3 , with Kza,o = (no preheating, 
dot-dashed), 0.4 (short-dashed), 0.8 (solid, best fits to the LT relation), and 1.6 (long-dashed). The observed distribution of 
the unidentified EGRET sources with |fe| > 45° is shown by the three solid lines, corresponding to all unidentified sources, 'em' 
sources (see Paper I), and steady unidentified sources defined by Gehrels et al. (2000) fll8| , with the order of the line thickness 
from thinner to thicker. The sensitivity limits of the EGRET and GLAST experiments are shown in the figure. The lower 
panel: the mean redshift, cluster mass (in units of 1O 16 M0), and angular radius (in degree) of gamma-ray clusters brighter 
than a given flux are shown. The line markings are the same as the upper panel. 
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